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^ ■ Abstract 

> : 

^ I The Gross-Neveu model with chemical potential is investigated as a low- 

energy effective theory of polyacetylene. In particular, we focus on the abrupt 
change in the features of electric conductivity such as sharp rise in the Pauli 
paramagnetism at dopant concentration of about 6%. We will try to explain it 
by the finite density phase transition in the Gross-Neveu model. The thermo- 
dynamic Bethe ansatz is combined with the large-N expansion to construct 
thermodynamics of the Gross-Neveu model. A first-order phase transition 
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is found in leading order in the 1/N expansion and it appears to be stable 
against the 1 /N correction. The next to leading order correction to the critical 
dopant concentration is explicitly calculated. 
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I. WHAT IS INTERESTING IN ELECTRIC CONDUCTIVITY IN 

POLYACETYLENE? 



In this talk we want to describe our recent works on thermodynamics of the Gross- Neveu 
model with chemical potential. The physical motivation for these works comes from 
condensed matter physics; electric conductivity properties of doped polyacetylene. 0] In 
particular, we are trying to describe the abrupt change in electric conductivity properties 
at dopant concentration of about 6%. 

Polyacetylene is a typical 1-dimensional polymer and we are particularly interested in 
trans-type polyacetylene. It is a fascinating material, showing variety of features of electric 
conductivity depending upon the dopant concentration. |^ At zero doping it is an insulator. 
Because of the Peierls instability an electron gas in 1 dimension is unstable against generation 
of charge density wave which entails, in this system, an alternating distortion of the 
lattice. It is called dimerization in condensed-matter literatures and the ground state is 
doubly degenerate. As a consequence an energy gap opens above the Fermi sea which 
implies an insulator. 

At dopant concentration up to ~6% the electric conductivity grows as dopant concen- 
tration increases but the Pauli paramagnetism stays low, indicating that the charge carrier 
is spinless. (A note for non-condensed matter physicists; doping is a procedure of importing 
the impurity atoms into the material to supply extra electrons or holes into the system. The 
dopant concentration is defined as the number of doped electrons per carbon atom.) It is 
believed that electric charge is transported by the spinless charged solitons in polyacetylene 
in this intermediate doping region. [|| One can say that the system displays the particular 
type of charge-spin separation. 

It appears that there is a consensus that the relevant dynamical degrees of freedom are 
solitons in this regime. Their existence is in fact predicted by Su, Schrieffer and Heeger 
in their original paper. 0] They introduced the SSH Hamiltonian to describe polyacetylene 
and discussed its physical implications. The spectrum of solitons in this model contains two 
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spinless charged solitons and two neutral solitons with spin 1/2. Prior to their pioneering 
work it was observed by Jackiw and Rebbi that in theories with fermions sohtons can 
carry fractional fermion quantum numbers. The peculiar features of spectrum of solitons in 
the SSH model can be nicely explained by the phenomenon of fermion fractionization. 

Then, at dopant concentration of about 6% an abrupt rise occurs in the Pauli para- 
magnetism while the conductivity keeps monotonically increasing. In some samples the 
conductivity reaches as high as that of copper. It is likely that the behavior indicates a 
transition to a metallic phase. This is the behavior that we want to understand. 

II. THE GROSS-NEVEU MODEL AS A LOW-ENERGY EFFECTIVE THEORY 

OF POLYACETYLENE 

Soon after the proposal of the SSH Hamiltonian, it was shown by Takayama, Lin-Liu and 
Maki PI that the low-energy effective continuum theory of the SSH model is given by a 1+1 
dimensional relativistic four-Fermi theory with two-flavor Dirac fermions. The Lagrangian 
of the model can be written by using an auxiliary fleld a (which represents the phonon 
degrees of freedom) and suppressing the flavor index as 

£ = ^27'^ - 1(7^ - gail)'^ (1) 

This is nothing but the fleld theory model which was flrst investigated by Gross and Neveu 
as an asymptotically-free renormalizable quantum fleld theory with many common features 
with quantum chromodynamics. 

We would like to utilize the Gross-Neveu model with chemical potential as the effective 
theory for a transition from solitonic to metallic phase. If it does a good job the model should 
display a phase transition at a certain value of the chemical potential. Is it likely? We believe 
that the answer is yes. It is believed that at zero chemical potential the ground state of the 
model violates the discrete chiral symmetry and it is doubly degenerate. It is, of course, in 
accord with the dimerized ground state and is nothing but the manifestation of the Peierls 
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instability. Therefore, the theory admits sohtons (kinks) and they are "fermion-fractionized" 
sohtons in the sense of Jackiw and Rebbi. The evidence for the chiral symmetry breaking 
and the doubly degenerate ground state comes from the large-N approximation. M It is also 



strongly supported by the existence of consistent factorizable fermion-fermion, ]T0|] fermion- 
kink and kink-kink S-matrix. | ]rT| , p!^ 

At large enough /i, on the other hand, we believe that the symmetry is restored. Because 
the Gross-Neveu model is an asymptotically free theory the fermions become non-interacting 
at high enough densities and are not able to support their condensate that they formed to 
break chiral symmetry at low densities. Therefore, it is very likely that the theory displays 
the phase transition at a certain fermion density or at the corresponding value of chemical 
potential. We note that this is in disagreement with what the soliton lattice theory predicts. 

HI 



III. THE GROSS-NEVEU MODEL WITH CHEMICAL POTENTIAL; THE 

LARGE-N LIMIT 

To verify (or disprove) this conjecture and to know the order of the phase transition 
(if any) we examined the Gross-Neveu model with chemical potential by using the large-N 
approximation. In the leading order of 1/N expansion the effective potential of the model 
can be obtained analytically by taking into account the effect of chemical potential. It reads: 



2N ^ AN 
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(a' - 7'){/n(aVao^) - 3} + ^(7' - a'){2ln ^ ^ - 3} 



2N 



fi^^i^'-a') (2) 



where 7 = /xy N/ vr. By analyzing the expression one can easily figure out that at jj = m/ v^, 
where m is the mass of the fundamental fermion, the theory has a triple degenerate ground 
states. (See Ref. 1 for details.) This is a clear signature for the first-order phase transition. 
One can evaluate the corresponding value of the fermion density, which can be interpreted 



as dopant concentration, as 

N 

where ^ is the correlation length (soliton size) measured in units of lattice constant. It is 
estimated as ^ ~ 7 by SSH. [|| It turns out that i/c = 0.064, very close to 6% at = 2. 
Here is the number of flavor of Dirac fermion in the effective field theory. The two flavor 
originates the spin degrees of freedom of the electrons in the SSH model. 

IV. GOING BEYOND THE LARGE-N LIMIT 

It is a good news that the large-N works but there arise a number of questions. 

(1) Are there any experimental evidences for the first-order phase transition? 

(2) Is the large-N approximation reliable even at N=2? 

(3) Can the result be viewed as a robust evidence for the presence of finite density phase 
transition in the Gross-Neveu model? 

Again a good news about (1) is that apparently there exists an experiment that signals a 
hysteresis in the Pauli paramagnetism as a function of chemical potential; a clear signature 
for first-order phase transition. The result, however, seems to be neither confirmed nor 
refuted by the other independent experiments. 

The second and the third questions are more pressing to us as theorists. Certainly we 
should be able to do better. As a first step toward answering these questions we have 
calculated the next to leading order correction to see if the first-order phase transition 
survives (we mean, if there are any instabilities occurring) and to compute correction to the 
critical dopant concentration. 

To carry this out, however, we had to invent a novel way of computation because the 
direct calculation of the effective potential at two-loop with chemical potential looks to be 
formidable. The method we employed is a hybrid method combining the thermodynamic 
Bethe ansatz (TEA) and the 1 /N computation at zero chemical potential. 
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V. THE THERMODYNAMIC BETHE ANSATZ 



The thermodynamic Bethe ansatz [^,0] is a powerful method for constructing thermo- 
dynamics of a 1 dimensional system of particles whose S'-matrix elements Sab are exactly 
calculable. The technique is particularly suited to our problem because all the S'-matrix 
elements of the Gross-Neveu model are known. UToHT^ By solving the integral equation of 
TBA one can, in principle, construct an exact thermodynamics of the Gross-Neveu model. 

Let us restrict ourselves to the elastic (or diagonal) scattering theories in which neither 



inelastic nor "charge-exchange" reactions occur. ||T^ Following Yang and Yang one can 
write down the integral equation for dressed energy ea{6) of quanta as a function of the 
rapidity 6, where the subscript a specifies the particle species. It reads, 

Ttaie) = -fia + ma cosh ^ - ^ / dO'KabiO - 9')ln [l + e-^^(^')] (3) 

b 

where the kernel Kab{0) is defined by 

Having solved the equation for taiO) the free energy density of the system can be computed 



as 



fil^) = -J2 [ dOrUa cosh ein [l + e"^"] + ^ ^aDa (5) 

where Da denotes the number density of particle a. We note that (negative of) the first term 
in (^) implies the pressure P because of the thermodynamic relation / = fiD — P. Hereafter, 
we consider the case of single species of particles. We argue that the approximation can be 
justified in the large-N expansion of the Gross-Neveu model. 

Since we are interested in the possibility of phase transitions at finite chemical potential 
we want to take zero temperature limit. It is because there is no phase transition at finite 
temperature in 1+1 dimensions due to the Mermin- Wagner theorem. |jl8| R One can show 



^ One can argue in length that the the Gross-Neveu model at zero-temperature serves as a better 
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that by taking zero-temperature limit T — with Te{6) = —e{6) kept finite, the TBA 



equation for "energy density" takes the form, |19 



e(^) = /X - mcoshe + r de'K{e - e')i{9'). (6) 

J — B 

where B is determined by e{B) = and the integration region is restricted to i{B) > 0. 
Using the solution of this integral equation the pressure can be expressed in zero-temperature 
limit as 



, m 

'2^' J-B 



^ d9coshee{e) (7) 



The negative of the pressure is the relevant quantity to discuss the finite-density phase 
transition and to calculate the critical dopant concentration. It is related with the grand 
canonical partition function S as —PL = —TlnS. ||2^ (Recall the similar relationship F = 
—TlnZ between the Helmholtz free energy F and the canonical partition function Z.) The 
fermion number density n is then given by — 

There are, however, obstacles both technical and conceptual. The technical problem is 
that the integral equation is difficult to solve analytically. But it can be done numerically. 
On the other hand the conceptual problem is more difficult to handle. We do not know how 
the phase transition can be signalled in the framework of the TBA. 



VI. THE TBA IN LEADING ORDER IN 1/N EXPANSION 

To gain insight to this problem we have examined the question by using the 1/N expan- 
sion. At large N one can argue that the only relevant degrees of freedom are fundamental 



effective theory of pure sample of polyacetylene at the room temperatures. The Mermin- Wagner 
theorem holds because of a strict one-dimensionality of space and it is known that its prohibition 
becomes invalid upon turning on very tiny three-dimensionality such as infinitesimal interchain 
couplings. H. M. thanks Yutaka Okabe for discussion on this point. 
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fermions because everybody else is more massive, with the masses being proportional to N. 
Therefore, the single species approximation should be valid. We have found an interesting 
structure. In leading order of 1/N expansion we have the analytic expression (0) of the effec- 
tive potential of the Gross-Neveu model with chemical potential. From this we can obtain 
the absolute value of the ground state energy density as a function of chemical potential: 

T/ /)/ 2 2^^^, + VfJ''^ - //^^ f^s 

Veff = — : V Oiu — m )-—\ln \ — ^ — 1. (8) 

Here we have chosen the renormalization point (Tq so that A = g'^N = tt and then ctq gives 
minimum of the effective potential. 

On the other hand, we can compute the free energy of the model by solving the TEA 
equation. In leading order of 1/N it is trivial to solve the equation because the kernel K{9) 
vanishes (namely, S = 1). We obtain for (negative of) the pressure, —P, the exactly the 
same expression as the second term in (^. The pressure given by the TEA is normalized 
such that it vanishes at = m and there is no way (to our knowledge) of computing it below 
= m. 

It is conceivable that the pressure — -P(/i) has additional /i-independent contribution 
— P(0). It can be interpreted as the vacuum energy density of the Gross-Neveu model. On 
dimensional ground it can be written as — -P(O) = —hiv?, where 6 is a constant. While 
there is no way of computing h within the framework of the TEA it is calculable, at least 
in principle, in the 1/N expansion. In the leading order it was computed by Gross and 
Neveu @ and h = l/An. Together with the TEA result of —P{fi), it reproduces the exact 
expression of the free energy, i.e., the effective potential (||) obtained by the explicit 1/N 
calculation. 

This establishes our strategy of how to compute free energy of the Gross-Neveu model 
order by order in 1/N expansion. Namely, we obtain —P{fi) by solving the TEA equation 
and supplement it by the vacuum energy density calculated at zero chemical potential. 

Then, the question is how the TEA free energy can signal the phase transition? We 
propose to take massless limit of the pressure —P^ji) — -P(O) calculated by the above method 
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to obtain the pressure in the massless phase. p[ We have 

ZTT 

the free fermion result. Then, we ask if the two pressures cross at certain value of the chemical 
potential. They do at /i = m/v^, reproducing the phase transition point predicted by the 
explicit leading large-N computation. 

VII. THE TBA IN NEXT TO LEADING ORDER 

In our second paper we have gone one step further by generalizing the procedure to the 
next to leading order in 1/N expansion. The TBA equation is no longer trivial to solve 
but fortunately it can still be solved in a closed form. We skip all the technicalities and just 
quote the result. 0,|T^ The 1/N correction to the free energy of the Gross-Neveu model can 
be given in a parametric form: 

2 

- P(^) = -j^[B^ + sinh'^B - BShi{2B)] (10) 

/i = mcoshB + ^ [mcoshB\Shi{2B) - 2B] - msinhB[Chi{2B) - ln{sinh2B) - 7]] (11) 

where Shi{x) and Chi{x) are the hyperbolic integral functions defined by 

, sinht 



Shi{x) = / di- 
Jo 



t ■ 

Chi{x) = 7 + /nx + / dt^^^^^ — -. (12) 
Jo t 

The uncalculable constant, the vacuum energy density at zero chemical potential, was 
indeed calculated long time ago by Schonfeld. PT| With a bit of numerical computation the 



correction to the constant b can be written as 

Then, one can go through the exactly the same procedure as the leading order. We still 
obtain the free fermion result by taking the massless limit but with chemical potential with 
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1/N correction. It is what we expect because of the asymptotic freedom of the Gross- Neveu 
model. (See Ref. 2 for details.) The resulting correction to the critical point can be written 

as 

m , 0.47, 

Therefore, it gives rise to about 20% correction to Hc at N=2. If it is expressed by critical 
dopant concentration i/c — 0.05, which is somewhat smaller than but still roughly agrees 
with the experimental value. 



VIII. CONCLUSION AND OUTLOOK 

In this talk we have discussed how far one can proceed toward the thermodynamics of 
the Gross-Neveu model at zero temperature and finite background fermion density. At least 
in the 1/N expansion we have found that there is persistent first-order phase transition at 
a certain value of the chemical potential. Whether it survives in an exact treatment at 
N=2 theory remains to be seen. However, we have argued that it is the case by relying 
on the following robust two features of the model: one is the chiral symmetry breaking at 
zero chemical potential as manifested in the kink spectrum which is known exactly by the 
method of factorizable S-matrix, and the other, the asymptotic freedom of the theory. 

Is it possible to go beyond the large-N expansion? It is a very nontrivial question. In 
principle the answer is yes. Since all the S-matrix elements are known it should be possible 
to formulate the complete thermodynamics of the Gross-Neveu model. A more difficult 
but challenging question is whether the phase transition can be signaled solely within the 
framework of the TEA. We are now trying to answer (at least some of) these questions. 
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